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^ ' ABSTRACT 

H \ We present the 90/im luminosity function of the Final Analysis of the European Large 

Area ISO Survey (ELAIS), extending the sample size of our previous analysis (paper 
IV) by about a factor of 4. Our sample extends to z = 1.1, ~ 50 times the comoving 
volume of paper IV, and 10 7 7 < h^ 2 L/L Q < 10 125 . From our optical spectroscopy 
campaigns of the northern ELAIS 90/im survey (7.4 deg 2 in total, to Sgo^m > 70mJy), 
we obtained redshifts for 61% of the sample (151 rcdshifts) to B < 21 identified at 
7^m, 15^m, 20cm or with bright (B < 18.5) optical identifications. The selection 
function is well-defined, permitting the construction of the 90/im luminosity function 
of the Final Analysis catalogue in the ELAIS northern fields, which is in excellent 
agreement with our Preliminary Analysis luminosity function in the ELAIS SI field 
from paper IV. The luminosity function is also in good agreement with the IRAS-based 
prediction of Serjeant & Harrison (2004), which if correct requires luminosity evolution 
of (1 + z ) 3 4±10 for consistency with the source counts. This implies an evolution in 
comoving volume averaged star formation rate at z ~ 1 consistent with that derived 
from rest-frame optical and ultraviolet surveys. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: star-burst - infrared: galaxies - submillimetre 



1 INTRODUCTION 

Recent infrared surveys have revealed the existence 
of strongly evolving populations of far-infrared-luminous 



galaxies. For example, the discovery of sub-mm luminous 
galaxies (Smail et al. 1997, Hughes et al. 1998) with median 
redshifts > 2 (Chapman et al. 2003) proves the existence 
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of a radical change in the demographics of star formation. 
Luminous infrared galaxies (~ 10 11 L@) contributed at least 
as much to the comoving star formation rate at z ~ 2 as 
all optical/ultraviolet identified sources, but are negligible 
contributors locally. Furthermore, the spectrum of the ex- 
tragalactic background (e.g. Puget et al. 1996, Fixen et al. 
1998, Lagache et al. 1999) implies that approximately half 
the luminous energy emitted by stars throughout the history 
of the Universe was absorbed and re-emitted by dust in the 
far-infrared. Differences between the evolution of infrared- 
luminous galaxies and that of galaxies selected in the rest- 
frame optical/UV may be in part attributable to the changes 
in the demographics of obscured vs. unobscured star forma- 
tion, but may also be related to large scale structure density 
variations in the small-area surveys (e.g. <C 1 square degree) 
conducted to date (see discussion in Oliver et al. 2000). Con- 
sequently, while some of the most significant earlier contri- 
butions have been in small-area surveys (e.g. Hughes et al. 
1998), the emphasis has since shifted to wide-area surveys 
(~ 1 — 100 square degrees) which are less affected by cosmic 
variance problems (e.g. Oliver et al. 2000, Lonsdale et al. 
2003). 

The European Large Area ISO Survey (ELAIS) was 
the largest open time survey on the Infrared Space Obser- 
vatory (ISO), and resulted in the delivery of the largest cat- 
alogue of any survey on ISO (Rowan- Robinson et al. 2004) 
from both the ISOCAM instrument (Cesarsky et al. 1996) 
and ISOPHOT instrument (Lemke et al. 1996). The science 
goals are diverse, and discussed in detail in Oliver et al. 
(2000, paper I). Following the Preliminary Analysis of the 
ELAIS data (Serjeant et al. 2000, paper II, Efstathiou et al. 
2000, paper III), a more exhaustive analysis of the data has 
been made resulting in deeper and better calibrated source 
lists (Lari et al. 2001, Gruppioni et al. 2002, Heraudeau et 
al. 2004 paper VIII), and is referred to as the ELAIS Final 
Analysis. The ELAIS survey complements the deeper sur- 
veys conducted by ISO (e.g. Franceschini et al. 2002 and 
references therein; Rodighiero et al. 2003) and the shallower 
ISOPHOT Serendipity Survey (Stickel et al. 2004). 

The Final Analysis catalogue paper (Rowan- Robinson 
et al. 2004) presents the source list and redshifts obtained to 
date from ELAIS. These sources have been studied in ear- 
lier papers in this series. For example, Morel et al. (2001, pa- 
per VI) present the first hyperluminous galaxy discovered in 
ELAIS, several more of which are catalogued in the Rowan- 
Robinson et al. (2004) paper; Pozzi et al. (2003) compare the 
15/im, Ha and 1.4GHz star formation rates of galaxies in the 
S2 field; Vasainen et al. (2002) report the near-infrared iden- 
tifications of ELAIS galaxies, in the process obtaining a im- 
proved flux calibration which implies a modest correction to 
the Preliminary Analysis source counts (Serjeant et al. 2000) 
consistent with the Final Analysis source counts (Gruppioni 
et al. 2002); Basilakos et al. (2003, paper VII); and Man- 
ners et al. (2003a, b) examine the X-ray properties of ELAIS 
galaxies. The evolution of AGN in ELAIS is examined in 
Matute et al. (2002). The ELAIS fields have also been the 
subject of many subsequent multi-wavelength surveys (e.g. 
Ciliegi et al. 1999, Gruppioni et al. 1999, Alexander et al. 
2001 (paper V), Scott et al. 2002, Fox et al. 2002, Manners 
et al. 2003a, b) making the fields among the best studied 
degree-scale areas on the sky. Finally, the main ELAIS sur- 
vey areas will shortly be observed by the SWIRE Legacy 



Survey (Lonsdale et al. 2003) on the Spitzer Space Tele- 
scope, making this catalogue extremely important for the 
rapid science exploitation of the Spitzer data. 

The multiwavelength galaxy luminosity functions are 
one of the main observational constraints on the evolution 
of dust-enshrouded star formation in galaxy populations. In 
Serjeant et al. (2002, paper IV) we derived the 90^tm lumi- 
nosity function of ELAIS galaxies in the Preliminary Analy- 
sis of the ELAIS SI field (Efstathiou et al. 2000). We found 
good agreement with local determinations at the low redshift 
end, and evidence for pure luminosity evolution at the rate of 
(1 + z) 2 ' 45±0,85 . This is consistent with the optically-derived 
evolution in the comoving volume-averaged star formation 
rate (e.g. Glazebrook et al. 2003). In this paper we present 
the 90/im luminosity function for the Final Analysis of the 
remaining areas (Heraudeau et al. 2004, paper VIII). Sec- 
tion [5] describes the sample selection and the data acquisi- 
tion, and section |3 presents the 1/Vmax luminosity function. 
SectionQJdiscusses our results in the context of source count 
model evolution and the comoving star formation history. 

We assume a Hubble constant of Ho = lOO/i = 72 
km s _1 Mpc -1 , and a cosmology of Qm = 0.3, Qa = 0.7 
throughout this paper. We adopt the convention of convert- 
ing from 90/xm monochromatic luminosities to bolometric 
luminosities assuming vL v =constant, i.e. 10 10 Lq corre- 
sponds to approximately 1.159 x 10 24 W Hz -1 at 90^im. 
This choice simplifies the comparisons with other studies 
and does not affect the results in this paper. 



2 METHOD 

2.1 Sample selection 

Our principal sample is the ELAIS 90/im Final Analysis 
sample of Heraudeau et al. (2004) in the Northern ELAIS 
fields, covering a total of 7.4 square degrees to a flux limit 
of 70mJy. The completeness as a function of flux is given in 
Heraudeau et al. (2004), and varies from ~ 100% at 150mJy 
to ~ 50% at 70mJy. The source counts at the faintest end 
(e.g. 70 — 130^iJy) differ from those of the Preliminary Anal- 
ysis (paper II, paper IV). These data points were neglected 
in the fit to the source counts, from which the evolution 
rate was derived, so the evolution quoted in paper IV is un- 
affected by the change in the catalogues at the faintest end. 
Notably, the source counts of the optically identified sources 
are identical between the Preliminary and Final analyses. 
We do not have Final Analysis catalogues for the small 
ELAIS survey areas (paper I), so we supplement the Fi- 
nal Analysis catalogue with optically identified sources from 
the Preliminary Analysis, which are almost entirely in these 
supplementary small fields. The exclusion of these sources 
does not affect determined luminosity function, apart from 
decreasing the spectroscopic sample size by seven objects. 
The catalogue has a total of 420 galaxies, of which we have 
optically identified 249 to magnitude limits discussed be- 
low. We have obtained optical spectroscopy of 151 of these 
optically identified galaxies to date. 

Figure plots the luminosity-redshift plane for our 
sample. For the purposes of K-corrections, the rest-frame 
galaxy spectrum in the vicinity of 90^im is assumed to fol- 
low dlog L„/dlog v = —2, though this choice has very little 
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Figure 1. Luminosity-rcdshift plane for spcctroscopically iden- 
tified ELAIS 90/im Final Analysis galaxies from the northern 
ELAIS fields. Monochromatic luminosities at 90^tm have been 
converted to bolomctric luminosities assuming i^L„=constant. No 
correction has been made for pure luminosity evolution. 

affect on the derived luminosity function. There are no pho- 
tometric redshifts used in the derivation of the luminosity 
function in this paper. 

2.2 Data acquisition 

ELAIS 90/im galaxies were optically identified either purely 
from the Digitised Sky Survey (for which paper IV found 
B < 18.5 is the limit for reliable identifications), or from 
the abundant multi-wavelength data from the ISOCAM 
mid-infrared ELAIS data (paper II, Gruppioni et al. 2003) 
and/or the VLA follow-up mapping (Ciliegi et al. 1999). 
For the 90^m galaxies identified with ISOCAM or the VLA, 
optical identifications brighter than B < 21 were sought us- 
ing the Digitised Sky Survey and our own optical imaging 
(Rowan- Robinson et al. 2004). Given the astrometric errors 
of the multiwavelength identifications (~ 1" — 3") this mag- 
nitude limit is much more than sufficiently bright to ensure 
reliable optical identifications. 

The optical identifications themselves are presented and 
discussed in Rowan- Robinson et al. (2004). It is possible that 
despite the careful analysis in that paper, some of the iden- 
tifications may be wrong; nevertheless, the almost-universal 
appearance of emission lines in our spectra (see below) im- 
plies that the overwhelming majority of the identifications 
must be correct. Even in cases where the identification is 
wrong, the spurious identification will often have the same 
redshift as the true identification since star-forming galax- 
ies well-known to be more likely to have close companions 
than the field galaxy population as a whole (e.g. Surace et 
al. 2004). Similarly, it is possible that at the ~ 40" reso- 
lution of ISO, some point sources will in fact be blends of 
multiple sources. However, the ELAIS 90/im survey is far 
from the point source confusion limit (e.g. Rowan-Robinson 



Figure 2. Completeness of the optical spectroscopy, as a func- 
tion of apparent B magnitude. At each magnitude m, the spec- 
troscopic completeness plotted refers to galaxies with magnitudes 
in the interval m — 1 tom + 1. This curve has structure on scales 
smaller than ±1 magnitude (e.g. discrete changes) as individual 
objects move in and out of the bin being considered. 

2001), so any blends are likely to be physically associated 
galaxy pairs, such as interacting or merging systems. Again, 
this would not lead to incorrect redshift estimates for the 
purposes of the luminosity function. Also, in local galaxies 
there is little evidence that galaxies in pairs each contribute 
comparably to their total unresolved far-infrared flux (e.g. 
Surace et al. 2004). In summary, for the purposes of this 
paper we will assume that the identifications are all correct. 

Spectroscopic campaigns were made on several runs at 
the William Herschel Telescope; the Guillermo-Haro 2.2m, 
Cananea, Mexico; and the Kitt Peak WIYN facility. The 
Kitt Peak and William Herschel data is described in Perez- 
Fournon et al. (in preparation). The Guillermo-Haro data 
was taken on 17-30 June 2001, 6-13 May 2002 and 1-10 July 
2002 using the LFOSC and Boiler & Chivens low dispersion 
spectrographs. Redshift errors are dominated by wavelength 
calibration accuracy (~ 1 — 2 A) and are negligible for the 
purposes of this paper. In total, we obtained spectroscopic 
redshifts of 61% of the sample, listed in Rowan-Robinson et 
al. (2004). 



3 RESULTS 

A luminosity function can be constructed for any sample 
provided that the selection function can be accurately stated, 
and that there are no populations of objects which are un- 
detectable at any redshift. For a single flux-limited sample, 
the number density in a given luminosity bin is given by 

* = SV^ (1) 
with an associated RMS error of 

A* = y/w^ (2) 
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Figure 3. Shaded area shows the ±1<t luminosity function con- 
straint at 90/mi from the ELAIS Final Analysis of the northern 
ELAIS fields. Pure luminosity evolution of (1 + z) 3 is assumed, 
though this has only a small effect on the derived luminosity func- 
tion except at the highest luminosities. Luminosities are converted 
to bolometric luminosities assuming vL v =constant. At each lu- 
minosity, the objects with luminosities ±0.5 dex are used in the 
construction of the luminosity function, which is indicated on the 
diagram as the boxcar width. Also plotted as the unshaded bound 
area is the 90/tm ELAIS luminosity function from the Preliminary 
Analysis of the SI field, which used the same boxcar width. Note 
the excellent agreement between the Preliminary and Final Anal- 
yses, and the much tighter constraint from the Final Analysis. 

where the sums are taken over the objects in the bin. The 
maximum volume Vm ax for each object is the comoving vol- 
ume enclosed by the maximum redshift at which such an 
object is detectable, given the flux limit. 

This procedure can easily be generalised for the pres- 
ence of more complicated selection criteria, including the 
multiwavelength flux limits used in this paper. For example, 
if the selection function is a function of redshift, then this 
function must be used to weight the differential volume el- 
ements, which are then integrated to calculate a weighted 
Vmax for each object. Paper IV (Serjeant et al. 2001) dis- 
cusses in more detail the formalism for generating 1/V max 
luminosity functions in the presence of complicated selection 
functions. In this paper we largely restrict our discussion 
to the construction of the selection function, and refer the 
reader to paper IV for its application to the luminosity func- 
tion. The only modification of this approach of Serjeant et al. 
(2001) to this paper is in the consideration of bins contain- 
ing only one galaxy. In this paper these bins (and only these 
bins) are treated as having a la range of (0.18 — 3.3)Vmax 
(rather than V^ax ± V^ x ), corresponding to the ±ler likeli- 
hood limits of an underlying Poisson distribution with one 
detection (Gehrels 1986). 

The optical completeness correction is the largest sys- 
tematic uncertainty in our selection function. To estimate 
this, we plot in figure[2]the fraction of galaxies with spectro- 
scopic redshifts, as a function of apparent magnitude. Note 



n ELAIS 90// in luminosity function 




L/L, 



Figure 4. Shaded area shows the ±ltr luminosity function at 
90/im from the ELAIS Final Analysis of the northern ELAIS 
fields, with a boxcar width of ±0.25dex (c.f. figurel^l. Pure lumi- 
nosity evolution of (1 + z) 3 is assumed. Luminosities are converted 
to bolometric luminosities assuming vL v =constant. Also plotted 
is the predicted 90/im local luminosity function from the PSC-z 
survey (Serjeant & Harrison 2004). Note the good agreement of 
the ELAIS luminosity function with this prediction. 



that this is a differential plot, i.e. at a magnitude m we only 
consider galaxies with magnitudes mil. We found that our 
luminosity function is robust to changes in the derivation of 
this optical completeness, such as considering galaxies with 
magnitudes m ± 1.5 or m ± 2 instead of m ± 1, in that any 
systematic changes to the derived luminosity function are 
very much smaller than the ±l<r random errors. 

To have an optical identification, each galaxy must ei- 
ther have a bright optical ID (B < 18.5) or be detected 
at 15/im, 7/im or in the radio. In the case of 15/im, 7/im 
and radio detections, the galaxy can be identified to fainter 
optical magnitudes (B < 21), since the greater astrometric 
accuracy of these catalogues makes reliable identifications 
possible to fainter magnitudes. The probability of a given 
galaxy having an identification is therefore given by 

Pit) = 1 - (1 - Popt)(l - Pl5nm)(l - P7 M m)(l - Prad) (3) 

Here p op t is the probability that this galaxy has a bright 
optical ID, and similarly for 15/im, 7/im and radio. For 
a spectroscopically-complete sample, p opt is 1 at redshifts 
where the galaxy would have B < 18.5, and at higher red- 
shifts. In our case, we use the optical completeness curve 
of figure for the B < 18.5 range. Similarly, pisjjm in- 
corporates both the 15/im completeness (e.g. Gruppioni et 
al. 2003) and the optical completeness from figure [5] i.e. 
Pisum = ci 5flm x c op t, where ci 5Mm is the 15/im complete- 
ness and c opt is the completeness of optically-faint galaxies 
in the optical. Analogous optical completeness corrections 
were made for 7/im and radio identifications. We calculate 
Pm(z) for every galaxy in our sample, and use it to weight 
the volume elements in the calculation of the total acces- 



ELAIS paper IX: the 90//m luminosity function 5 



90/im differential source counts 



z;, 

ID 



\ 

a 



10.0 P 



1.0 



0.1 




0.01 0.10 1.00 10.00 100.00 
90/i,m flux (Jy) 



Figure 5. Source counts from ELAIS (filled circles, Heraudeau 
et al. 2004), the Lockman Hole (open circles Rodighiero et al. 
2003), and local counts estimated from IRAS (crosses, Serjeant & 
Harrison 2004). Also overplotted are three model pure evolution 
curves of (1 + z) a , where a = (bottom curve, no evolution), 
a = 3 (central curve) and a = 4 (upper curve). All evolution is 
assumed to have stopped in these models at z = 2. Note the clear 
deviation of the source counts from the no-evolution case. 



Luminosity bin 



<1> 

Mpc" 



dcx 



Mpc" 



dcx 



8.625 ± 0.25 


1 


86 x 10 


_3 


1 


84 x 10" 


3 


8.75 ± 0.25 


1 


1 x 10" 


3 


2 


62 x 10" 


3 


8.875 ± 0.25 


8 


4 x 10" 


3 


3 


64 x 10" 


3 


9. ± 0.25 


1 


14 x 10 


_2 


i 


13 x 10" 


3 


9.125 ± 0.25 


1 


1 x 10" 


2 


3 


89 x 10" 


3 


9.25 ± 0.25 


9 


38 x 10 


_3 


3 


41 x 10" 


3 


9.375 ± 0.25 


6 


3 x 10" 


-i 


2 


35 x 10" 


3 


9.5 ± 0.25 


1 


07 




10 


_3 


1 


3 x 10" 3 




9.625 ± 0.25 


3 


17 


X 


10 


_3 


7 


84 x 10" 


4 


9.75 ± 0.25 


3 


88 


x 


10 


_3 


7 


68 x 10" 


4 


9.875 ± 0.25 


3 


02 


x 


10 


_3 


5 


57 x 10" 


4 


10. ± 0.25 


2 


75 


X 


10 


_3 


i 


73 x 10" 


4 


10.125 ± 0.25 


2 


06 


x 


10 


_3 


3 


49 x 10" 


4 


1 n or I a or 

10.25 ± 0.25 


1 


53 


x 


10 


_3 


2 


38 x 10" 


4 


1 375 + 25 


1 


09 


x 


10 


-3 


1 


87 x 10" 


4 


10.5 ±0.25 


6 


66 


x 


10 


-4 


1 


24 x 10" 


4 


10.625 ± 0.25 


3 


16 


x 


10 


-4 


6 


8 x 10" 5 




10.75 ±0.25 


1 


41 


x 


10 


-4 


3 


45 x 10" 


5 


10.875 ± 0.25 


5 


85 


x 


10 


-5 


1 


66 x 10" 


5 


11. ± 0.25 


2 


11 


x 


10 


-5 


9 


15 x 10" 


6 


11.125 ± 0.25 


1 


38 


x 


10 


-5 


7 


27 x 10" 


6 


11.25 ±0.25 


1 


38 


x 


10 


-5 


7 


4 x 10" 6 




11.375 ± 0.25 


1 


28 


x 


10 


-5 


7 


39 x 10" 


6 


11.5 ±0.25 


3 


98 


x 


10 


-6 


2 


71 x 10" 


6 


11.625 ± 0.25 


2 


92 


X 


10 


-6 


2 


61 x 10" 


6 


11.75 ±0.25 


1 


54 


X 


10 


-7 


1 


53 x 10- 


7 


11.875 ± 0.25 


3 


34 


X 


10 


-8 


+ 


7.68 lf) 
2.74 A lu 





Table 1. Tabulated luminosity function with boxcar width ±0.25 
dex, as plotted in figure HI 



sible volume Vm ax (Avni & Bahcall 1983) for each galaxy, 
following the procedure in paper IV. 

What is the best way of representing the luminosity 
function? One can make parametric fits, but these are sub- 
ject to the assumptions in the parameterisation; alterna- 
tively, one can directly bin the luminosity function, but this 
loses information. Dunlop & Peacock (1990) used a 'free- 
form' fitting technique in their analysis of the luminosity 
function of radio-loud active galaxies. In this methodology, 
a selection of high-order polynomials are fitted to the data. 
The arbitrariness and variation of the fitting functions partly 
avoids the issue of assumptions in the functional form. Di- 
vergence in the fits indicates lack of constraint, but unfor- 
tunately convergence does not necessarily imply constraint. 
Also, the range spanned by the models does not correspond 
to e.g. a ±lcr constraint. In this paper we adopt an alterna- 
tive approach, similar to that used in Serjeant et al. (2001): 
instead of binning the data into discrete bins, we calculate 
the luminosity function at every luminosity L in a bin of 
logarithmic width ±<5 log 10 L. This has the advantage of re- 
turning a ±1(7 constraint at any luminosity, while explicitly 
controlling the balance between resolution (in luminosity) 
and loss of information with the choice of the boxcar size. 

Figure |3] shows the 90/im 1/Vmax luminosity function 
derived in this way. We compute the luminosity function in 
a logarithmic bin of width ±0.5 dex. Figure [3] also gives a 
comparison with the ELAIS Preliminary Analysis luminos- 
ity function from the ELAIS SI field (Serjeant et al. 2001), 
and the agreement is excellent. However, a disadvantage of 
such a wide binning of the data is that there are effectively 
very few independent data points (e.g. 4 for a data set span- 
ning four decades in luminosity) . This was necessary for the 



small spectroscopic sample (37) in paper IV, but the current 
sample of 151 is large enough to benefit from finer binning. 
Therefore in figure 2] we plot the 90/im luminosity function 
for a ±0.25 dex binning of the data. This figure also plots a 
comparision with the predicted 90/im local luminosity func- 
tion from the PSC-z survey (Serjeant & Harrison 2004), es- 
timated from the IRAS colours of each of the 15411 PSC-z 
galaxies. In both figure |3] and figure 01 the population is as- 
sumed to evolve with pure luminosity evolution at a rate 
of (1 ± z) 3 . This assumption does not affect the luminosity 
function determined from the data except at the brightest 
end (figure 0. The luminosity function is tabulated in table 

m 



4 DISCUSSION AND CONCLUSIONS 

Using differential volume elements weighted to the selec- 
tion function (section [3J| , the (V/Vmax) distributions hint 
at the existence of evolution in this population, but only 
at the ~ 0.6a level: assuming a no-evolution model yields 
(V/Vmax) = 0.523 ± 0.029, while assuming (1 ± z) 3 pure lu- 
minosity evolution yields (V/V ms , K ) = 0.482 ± 0.029. With 
this test, deviation from 1/2 is evidence for evolution (Avni 
& Bahcall 1983). However, this statistic is not the most ef- 
ficient at detecting evolution. 

Figure 2] demonstrates that the ELAIS 90/im luminos- 
ity function is consistent with the prediction from the lo- 
cal IRAS population by Serjeant & Harrison (2004). If we 
assume that this prediction is correct, then we can use 
the source counts to constrain the strength of the evolu- 
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tion. Figure plots the ELAIS 90^tm Final Analysis source 
counts (Heraudeau et al. 2004), together with the interpo- 
lated 90/im IRAS counts from Serjeant & Harrison (2004). 
We also plot the predicted source counts for the predicted 
90^im luminosity function with (l + z) a pure luminosity evo- 
lution to z — 2, with a — 0,3, 4. Clearly, the source counts 
are a strong discriminant of the strength of the evolution. 
Provided the local luminosity function is of this form, the 
constraint on the evolution is 2.4 < a < 4.4 at 68% con- 
fidence, or 1.3 < a < 4.6 at 95% confidence. Is z — 2 the 
most appropriate redshift to halt the evolution? There is ev- 
idence that the evolution in the populations at 15/im stops 
at z ~ 1 (Franceschini et al. 2001, Chary & Elbaz 2001, Xu 
et al. 2003) . Decreasing the maximum redshift to which pure 
luminosity evolution extends from z — 2 to z — 1 slightly 
alters our constraints to 2.4 < a < 4.6 at 68% confidence, 
or 1.3 < a < 4.8 at 95% confidence. This evolution is con- 
sistent with that derived by us in paper IV. Note that pure 
density evolution (PDE) can already be excluded as it over- 
predicts the sub-mjy radio source counts (Rowan- Robinson 
et al. 1993), and we do not consider PDE models further in 
this paper; however, it is worth noting that the far-infrared 
data alone is not sufficient to rule out pure density evolution. 

This evolution rate is in good agreement with that de- 
rived from optically-selected galaxy samples (e.g. SLOAN, 
Glazebrook et al. 2003). Our data is consistent with the 
obscured star formation history and the corresponding un- 
obscured history both evolving at similar rates, although 
the data at 90pm also permits stronger evolution than that 
determined from the optically-selected samples. This is in 
agreement with the determinations from the sub-mjy radio 
population (e.g. Haarsma et al. 2000, Gruppioni et al. 2001) 
which are sensitive to both the obscured and unobscured 
high-mass star formation. However, the apparent plateau in 
the cosmic star formation history above z — 2 indicated by 
both optically-selected galaxies (e.g. Thompson 2003) and 
in the sub-mm (e.g. Hughes et al. 1998, Scott et al. 2002, 
Fox et al. 2002) is difficult to reconcile with the current co- 
moving mass density in stars, Q, (e.g. Blain et al. 1999, 
Serjeant & Harrison 2004). For the infrared- luminous pop- 
ulation, this may be in part due to contamination of the 
far-infrared flux by cirrus heated by the interstellar radia- 
tion field (Efstathiou & Rowan- Robinson 2003, Kaviani et 
al. 2003). Alternatively, the initial mass function at high 
redshift may be skewed to high masses (e.g. Larson 1998, 
Franceschini et al. 2001, Serjeant & Harrison 2004). 

The Spitzer satellite has already launched successfully, 
and the prospects for improving on this analysis in the near 
future are excellent. The SWIRE survey (Lonsdale et al. 
2003, 2004) will survey 70 square degrees including all the 
major ELAIS fields, deeper at 70/im than the ELAIS survey 
at 90/^m. The source counts will be a powerful discriminant 
of the evolution, but more importantly the existence of large 
spectroscopic redshift surveys in the SWIRE fields (Rowan- 
Robinson et al. 2004) will be a powerful and immediate tool 
for the exploitation of this important legacy data set. 
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